Science China Life Sciences

Virus-like particles (VLPs) enable transient, non-integrating delivery of CRISPR-Cas9 ribonucleoprotein cargo. Although VLPs have been reported for efficient DNA editing via base editors RNP delivery, the diversity of base editors tested as VLPs remains limited. We generated and benchmarked a panel of 12 base editors on the v5 eVLP backbone, targeting three genomic loci (HEK3, B2M, PDCD1) across five VLP dosages in LentiX-293T cells. Editing efficiency was generally dosage-dependent across all editors and varied by editor class and identity; PAM-flexible variants had lower editing efficiency than NGG-restricted counterparts, and the dual-function SPACE base editors showed reduced efficiency. We further characterized position-specific editing efficiencies and outcomes of the base editor VLP collection, revealing that a wide variety of mutation types are possible with the base editors in this collection.

Although Cinnamomum includes globally important spice and medicinal plants, genome-scale knowledge of Cinnamomum species native to Japan remains limited. The Nansei Islands are an approximately 1,200-km-long archipelago located between the Pacific Ocean and the East China Sea and represent a key biogeographic region shaped by deep-sea barriers such as the Tokara and Kerama Gaps. Here, we integrated chloroplast and nuclear genome data generated in this study with public datasets to examine the phylogenetic relationships of four Cinnamomum species native to the Nansei Islands and ten related species. Chloroplast genomes showed no major structural differences among the four island species, whereas nuclear genome data resolved two distinct evolutionary lineages in the Nansei Islands: C. yabunikkei, C. x durifruticeticola, and C. daphnoides formed a distinct Nansei Islands lineage separate from Sri Lankan, Indonesian, and Taiwanese lineages, whereas C. sieboldii was placed in a Taiwanese lineage. Ancestry estimation indicated that C. x durifruticeticola and C. daphnoides share a common ancestral component, while C. sieboldii exhibited a single homogeneous ancestral component, suggesting reduced ancestral diversity after divergence from a Taiwanese lineage. These findings reveal previously unrecognized genomic differentiation among Japanese Cinnamomum species and suggest that deep-sea barriers, dispersal, and long-term isolation shaped their diversification.

Kolmioviruses (KoVs; family Kolmioviridae) are circular, single-stranded RNA viruses that have recently been discovered in a wide range of animal species. However, invertebrate-associated KoVs have so far been identified only in termites and their diversity and evolution remain largely unexplored. Here, we comprehensively analyzed publicly available animal RNA-seq datasets and identified six novel KoVs associated with one fish and five invertebrate species: coho salmon (Oncorhynchus kisutch [Walbaum, 1792]), spotted babylon (Babylonia areolata [Link, 1807]), Atlantic pearl oyster (Pinctada imbricata [Roding, 1798]), small emerald moth (Hemistola chrysoprasaria [Esper, 1795]), codling moth (Cydia pomonella [Linnaeus, 1758]), and cocoa pod borer (Conopomorpha cramerella [Snellen, 1904]).

Porcine deltacoronavirus (PDCoV) is a globally circulating coronavirus of swine that causes severe enteric disease and mortality in neonatal piglets. Its broad host range and zoonotic potential raises concerns about continued spillover into human populations. However, the host factors that govern PDCoV infection remain largely unexplored. Identifying host factors involved in virus infection may uncover novel therapeutic targets for antiviral drug development or targets for gene editing to create resistant pigs. Here, we conducted a genome-wide CRISPR knockout (GeCKO) screen to identify host factors contributing to PDCoV infection in swine testicular (ST) cells. Small guide RNAs (sgRNAs) targeting the ANPEP gene, expressing the main receptor porcine aminopeptidase N for PDCoV, were enriched in PDCoV-resistant cells, demonstrating the validity of the screen. Moreover, we identified transmembrane protein 30A (TMEM30A), an accessory {beta}-subunit of the P4-ATPase flippase complex, as a host dependency factor for PDCoV infection along with several potential PDCoV-induced cell death regulators, including DET1, IFNAR1, RFWD2, STAT1 and TP53. We further demonstrated that PDCoV attachment and internalization into cells were significantly impaired in TMEM30A knockout ST cells, suggesting that TMEM30A is important in the early stages of the PDCoV life cycle. Because TMEM30A also contributes to SARS-CoV-2 replication, future studies are warranted to investigate its mechanistic role in the context of coronavirus infection.

Background: Large-scale biobanks have advanced genomic and epidemiologic research, but many rely on infrequent biological sampling and limited digital phenotyping. The Yamanashi Multi-omics Cohort (YMoC) was established to support longitudinal assessment of molecular, clinical, and behavioural changes in a screening-defined cohort of adults at elevated metabolic risk without diagnosed diabetes. Methods: YMoC is a longitudinal cohort of 215 adults aged 30-70 years in Yamanashi Prefecture, Japan, who met prespecified glycaemic eligibility criteria at health check-up, including fasting plasma glucose 100-125 mg/dL (5.6-6.9 mmol/L) and HbA1c <6.5%. Participants underwent three in-person visits over six months. Measurements include 75-g oral glucose tolerance testing with serial sampling, clinical biochemistry, anthropometry, liver elastography, and collection of blood, urine, stool, and saliva for multi-omics profiling. Between visits, participants wore a Fitbit Inspire 3 and completed daily app-based questionnaires using the Taohealth app. Current molecular data include genome-wide single nucleotide polymorphism array genotyping and longitudinal plasma proteomics in a subset. Conclusions: YMoC is designed to evaluate within-person molecular and phenotypic trajectories in a screening-defined metabolic-risk cohort. The cohort provides a dense longitudinal resource linking clinical assessments, biospecimens, omics assays, and digital phenotyping, including analyses of insulin-resistance-related markers such as homeostasis model assessment of insulin resistance (HOMA-IR).

Begomoviruses present severe threats to global crop production through complex vector-mediated transmission by the whitefly Bemisia tabaci to host plants such as tomato (Solanum lycopersicum). Unraveling the molecular dialogue between host immune activation and vector non-coding RNA networks is essential for identifying key drivers of virus persistence and transmission. Public host transcriptomic (GSE309527) and vector small RNA (sRNA) sequencing datasets (GSE111343) were processed through a multi-omics harmonization and signal calibration pipeline. Differential expression analysis was performed using empirical Bayes moderated linear models, followed by non-parametric Spearman rank correlation modeling ({rho}) to infer cross-kingdom co-expression dynamics and pathway enrichment profiling across host and vector bio-systems. Harmonized principal component analysis showed clear separation by infection status across host plant and vector cohorts. Differential expression analysis identified 138 significantly altered host genes (69 upregulated, 69 downregulated) and 130 differentially expressed vector sRNAs (65 upregulated, 65 downregulated). Host responses were dominated by significant upregulation of gene-silencing machinery, including Suppressor of Gene Silencing 3 (SGS3; log2 FC = 3.67, q = 7.47 x 10-5), and pathway enrichment in Jasmonate defense (q = 0.0004) and RNA Interference & Silencing (q = 0.0001). Vector sRNAs exhibited targeted dynamic alterations, with pathway enrichment in Salivary Gland Secretion ($q = 0.0030) and Gut Endosymbiont Response (q = 0.0210). Cross-kingdom correlation modeling revealed two distinct, highly anticorrelated regulatory modules (mean |{rho}| = 0.76). Host SGS3 expression strongly correlated with vector sRNA VEC_0080 ({rho} = 0.9762) and virus-derived siRNA Bt-vsiRNA-01 ({rho} = 0.7619). These findings demonstrate a tightly synchronized tripartite molecular crosstalk between host antiviral immunity, viral siRNA accumulation, and vector small RNA remodeling. These cross-kingdom regulatory modules highlight promising targets for dual-action RNA interference strategies aimed at controlling Begomovirus transmission.

Compact Cas nucleases offer advantages over the widely used SpCas9 due to their smaller size, which enables more efficient delivery for in vivo applications. Among these, the phage-encoded Cas{Phi}2 (Cas12j2) is highly promising due to its relaxed PAM requirement (5-TTN-3) and compact size (757 aa); however, its translational potential is limited by low editing activity. To enhance the efficacy of Cas{Phi}2, we optimized the previously reported EPICA system, developing EPICA.2, a eukaryotic directed evolution platform to improve nucleases with nearly undetectable activity. EPICA.2 integrates additional yeast evolution rounds to enrich for active variants along with a low background mammalian reporter system that improves detection and selection of enhanced variants. Finally, we set up a long-read sequencing protocol which uses unique molecular identifiers (UMIs) to reduce sequencing errors, enabling accurate identification of the mutation combinations in each evolved variant. Among the most frequent variants, we obtained evoCas{Phi}2, which contains six activity-boosting mutations with a synergistic effect not predictable by rational engineering. Overall, evoCas{Phi}2 showed up to 70-fold increased activity in human cells compared to wild-type and outperformed variants generated through rational approaches, highlighting the potential of EPICA.2 as a powerful strategy to evolve genome editing tools with low native activity.

Giant viruses (GV) are increasingly recognized as important ecosystem regulators. While metagenomics has uncovered extensive GV diversity, the global distributions of individual species and the biogeographic processes driving the pattern remain poorly understood. Here, we reconstructed GV metagenome-assembled genomes (MAGs) from 35 globally distributed deep freshwater lakes spanning five continents, aiming to identify their biogeographic patterns. The resulting 1663 non-redundant MAGs significantly expanded the known freshwater GV diversity, with [~]84% lacking a previously reported species representative. These MAGs were grouped into cosmopolitan and geographically restricted lineages. We identified 27 cosmopolitan GV species spanning multiple viral lineages, including families of Imitervirales, Pimascovirales, and mirusviruses order Styxvirales. The cosmopolitan species were characterized by their larger genomes and expanded gene repertoires of host-interaction functions, which may facilitate interactions with diverse hosts and contribute to their global distributions. The presence of geographically restricted species and the stronger distance-decay in community similarity observed in freshwater than marine ecosystems suggest that physical connectivity between ecosystems is an important factor influencing GV dispersal. We identified 312 and 177 GV MAGs almost exclusively associated with the epilimnion and hypolimnion, respectively. This water-layer preference of individual MAGs was highly consistent across lakes, suggesting conserved vertical partitioning in association with the thermal stratification of the water column. Overall, our findings reveal that GV biogeography in deep freshwater lakes is structured by the combined influence of horizontal dispersal limitation, vertical partitioning, and lineage-specific evolutionary histories.

Parasite communities in the internal organs of small mammals remain poorly characterized because metatranscriptomic pathogen discovery has largely focused on viruses and bacteria. We individually sequenced liver transcriptomes from 351 small mammals collected across 14 cities in Guangdong Province, China, and combined sequence-based detection with phylogenetic, abundance, transcriptional, and morphological analyses. Parasites were detected in 148 animals (42.2%). Four liver-associated taxa were identified, including three putatively novel or highly divergent lineages. Trichinellidae sp1 dominated the community, occurring in 36.2% of animals (127/351), five host species spanning Rodentia and Eulipotyphla, and all surveyed cities. Hydatigera taeniaeformis occurred in 3.7% (13/351), and its identification was supported by hepatic cysts and morphologically recovered strobilocerci. Eucoccidiorida sp1 and Trypanosoma sp1 were detected sporadically. Mapping to available reference genomes revealed transcripts from multiple genes and functional categories, supporting genuine biological signals rather than isolated marker-gene matches. Amsterdam virus covaried with Trichinellidae sp1 (r = 0.59), whereas Wenzhou tapeworm virus 1 covaried with H. taeniaeformis (r = 0.47). These findings reveal an unexpectedly prevalent, phylogenetically diverse, and broadly distributed liver-associated parasite community and establish individual-organ metatranscriptomics as a useful framework for parasite discovery and surveillance.

A protein whose removal or deletion causes significant disruption or collapse of a protein-protein interaction (PPI) network is referred to as a vulnerable protein. Such proteins may serve as valuable therapeutic or diagnostic targets in disease-associated networks. In this study, two PPI networks were constructed, one for HPV-positive and the other for HPV-negative head and neck squamous cell carcinoma (HNSCC), and the vulnerable proteins of these networks were identified by the node deletion approach. After analyzing the networks, 27 unique vulnerable proteins in HPV-positive and 72 unique vulnerable proteins in HPV-negative HNSCC were identified. Among them, one HPV-positive and seven HPV-negative HNSCC vulnerable proteins were further chosen by integrating multi-omics data. To exploit the vulnerabilities of these proteins, candidate synthetic lethal (SL) partners were predicted whose inhibition may selectively impair tumor survival. Subsequently, drug-gene interaction analysis was performed to identify inhibitors targeting the SL partners of these vulnerable proteins. Notably, in HPV-positive HNSCC, TOP2A, CHEK1, and CHEK2 genes were identified as SL partners of TTN, and their inhibitors were already clinically approved. While in HPV-negative HNSCC, ADA and MMP19 were identified as an SL partner of LMO7; TMEM45B, CDH3, and ELF3 genes were identified as an SL partner of CGN; and ZNF433 was identified as an SL partner of FLNC. However, MMP19, ZNF433, and TMEM45B inhibitors were not reported. Thus, these vulnerable proteins, including their SL partners, provide novel avenues to explore and develop more efficient and precise therapeutic and diagnostic strategies.

A central question in developmental biology is how spatiotemporal embryo morphological events and cell differentiation are precisely coordinated to eventually develop into a mature organism. Previous research has shown cell positions and asymmetric division stages are related to cell lineage specification since morula stage. However, the detailed spatiotemporal dynamics and histories of embryonic cells and their relations with cell lineage specification remain incompletely understood from the onset of embryo development onward. My study on in-vitro embryos showed that continuous live cell tracking mapped community-like patterns in cell origin and differentiation. Cell positions, cell temporal factors (such as cell division order and cell cycle) and their histories intricately interacted over pre- and peri-implantation stages. These interactions of spatiotemporal cellular activities directed embryo morphological events and asymmetric cell lineage origin and differentiation throughout 2-cell to around 100-cell stage. In conclusion, this study provides comprehensive insights into the cellular spatiotemporal dynamics of cell inheritance and differentiation during pre- and peri-implantation morphological events and cell lineage specification. By longitudinally integrating spatial and temporal parameters relating to early embryogenesis, the findings not only bridge divergent explanations stemmed from different studies in the field but also offers a refined foundation for evaluating embryo potential and improving outcomes in assisted reproductive technologies, stem cell research, and regenerative medicine.

Desertification exerts dramatic selection pressures on the evolution of plants. Despite the key role of ecological adaptation by natural selection to arid grasslands and subsequent intraspecific divergence, specific mechanisms driving this process remain poorly understood. Psammochloa villosa, a perennial forage grass endemic to the arid grasslands in Northwest China, where it thrives in shifting and semi-fixed sand land due to its exceptional drought tolerance, provides an ideal system to study adaptive evolution to aridity. In our study, we assembled a high-quality, chromosome-scale genome and conducted genomic resequencing of 42 populations across its major distribution. The genome assembly, which is approximately 1.55 Gb in size, has a super-scaffold N50 of 66.79 Mb, with 75.84% of the sequences identified as transposable elements. Coalescent phylogeny and genomic collinearity analyses strongly supported that P. villosa and Neotrinia splendens, as the closest taxa, shared a recent whole-genome duplication (WGD) event occurring approximately 18-20 Mya and followed by their divergence around ~11.2 Mya. Based on ancestral grass karyotype (AGK) reconstruction from synteny analysis, our results suggest that, relative to the AGK after the {rho}-WGD event, P. villosa and N. splendens underwent similar chromosomal restructuring and lineage-specific retention of numerous copies, providing a genomic basis for potential ecological adaptation and intraspecific diversification. The expanded XTH family, which encodes enzymes mediating xyloglucan endotransglucosylation and hydrolysis and thereby regulating xyloglucan remodeling, showed strong transcriptional responses under PEG-6000 treatment, suggesting that retained copies may be associated with xerophytic adaptation in P. villosa. Together, these findings suggest that WGD-derived gene retention created a delayed reservoir of genetic diversity that was later shaped by desertification and Qinghai-Xizang Plateau environmental changes, contributing to climate-associated genomic islands and intraspecific differentiation.

Replication-competent coronaviruses carrying fluorescent protein-tagged structural proteins remain scarce. Using the highly attenuated pangolin coronavirus GX_P2V(short_3UTR) as a backbone, we generated GX_P2V-mBJ-N, a recombinant coronavirus in which the bright green fluorescent protein mBaoJin is fused to the nucleocapsid (N) protein. The reporter virus is attenuated and genetically unstable in normal Vero cells but can be amplified to high titers in cells expressing wild-type N, and its fluorescence directly reports N protein expression. Using this authentic-virus platform, we show that high-titer GX_P2V cross-neutralizing antibodies persist in most healthy individuals and that cepharanthine potently blocks viral entry. GX_P2V-mBJ-N thus provides a simple and reliable tool for coronavirus tracing, immune surveillance, and antiviral drug evaluation.

Cyclic oligosaccharides possess industrial advantages, including molecular encapsulation capability and high physicochemical stability, owing to the absence of a reducing end. Recently, a novel cyclic tetrasaccharide, cycloisomaltotetraose (CI4), consisting of four -1,6-linked glucose units, and the enzymes responsible for its synthesis, cycloisomaltotetraose glucanotransferases (CI4Tases), were discovered. Unlike known cycloisomaltooligosaccharide glucanotransferases (CITases) that yield a wide distribution of cyclic products with a degree of polymerization (DP) of 7 or higher, CI4Tases strictly produce CI4. To elucidate the molecular mechanism underlying this strict DP4 specificity, we determined the crystal structures of CI4Tase from Agreia sp. D1110, in its ligand-free form, as well as in complex with the linear hydrolysis product isomaltotetraose (IG4) and with CI4. Structural comparisons revealed that a loop (M247 to R251) blocks the region corresponding to the -5 subsite of typical CITases, narrowing the substrate-binding pocket. This "molecular ruler" mechanism ensures that only a glycan chain of exactly four glucose units is accommodated for cyclization. Among mutants of the residue positioned at the center of bound CI4, the formation of by-products other than CI4 was significantly suppressed in F245L, F245A, and F245W. While the cyclization activity of all F245 mutants decreased, the CI4 hydrolysis activity of these three mutants was also significantly reduced, resulting in an increased specificity for cyclic sugar production. These findings elucidate the strict size-control mechanism of CI4Tase and provide a structural foundation for engineering cycloisomaltooligosaccharide-producing enzymes with optimized transglycosylation efficiency and specificity for industrial applications.

Very short RNAs play critical roles in modern biology, and are thought to have been crucial for genome replication during the origin of life. Next-generation sequencing is an essential tool for characterizing pools of small RNAs, but current library preparation methods suffer from strong size and sequence biases. Here we present tinyRNA-seq, an optimized library preparation method designed to minimize length- and sequence-dependent capture bias enabling the sequencing of RNA fragments as short as 2 nucleotides. We use degenerate adaptor regions to reduce ligation sequence bias and facilitate unique molecular identifier (UMI) installation. We benchmarked tinyRNA-seq against commercial kits using a model primordial RNA genome consisting of hundreds of defined oligonucleotides ranging from 2 to 12 nucleotides. tinyRNA-seq reproduced the input RNA distribution without the size and sequence bias of the commercial kits. tinyRNA-seq also enables the detection of de novo oligonucleotide generation, an important process for the origins of life. Applied to biologically derived small RNAs including miRNAs, piRNAs, and cityRNAs, tinyRNA-seq showed significantly lower capture bias and recovered a wider range of sequences than commercial kits. tinyRNA-seq may thus provide a more complete and quantitatively accurate representation of small RNAs from both biological and chemical sources. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=94 SRC="FIGDIR/small/743385v1_ufig1.gif" ALT="Figure 1"> View larger version (29K): org.highwire.dtl.DTLVardef@95ee64org.highwire.dtl.DTLVardef@155fb06org.highwire.dtl.DTLVardef@1d3665forg.highwire.dtl.DTLVardef@1e61404_HPS_FORMAT_FIGEXP M_FIG C_FIG

The long-range RNA-RNA interaction between the termini of a flavivirus genome is controlled by multiple local RNA secondary and tertiary structures. However, a complete picture of this essential genome cyclization process is still elusive. Herein, by using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE), SHAPE and mutational profiling, in vitro RNA binding assays, and reverse genetics approaches, we demonstrated that the hairpin 2 in the capsid-coding region (cHP2) is required for efficient replication of yellow fever virus and functionally interacts with the 3'-small hairpin (sHP) during genome cyclization. Interestingly, through a combination of multidimensional chemical probing and functional characterization assays, it was discovered that the sHP forms previously unknown tertiary interactions with an immediately upstream GACG region embedded within the 3'-cyclization sequence. This crucial tertiary RNA motif accounts for a potential negative regulatory mechanism of genome cyclization and is universally conserved among the mosquito-borne-related flaviviruses, which constitute a major phylogenetic branch within the subgenus Euflavivirus, genus Orthoflavivirus. Our study highlights the intricacy of RNA-based regulation in flaviviruses.

Epstein-Barr virus (EBV) infects 90% of the population worldwide and is associated with multiple malignancies, including Burkitt, Hodgkin and post-transplant lymphomas, nasopharyngeal and gastric carcinoma. EBV also triggers autoimmune diseases, including multiple sclerosis and systemic lupus erythematosus. While host factors important for EBV B-cell entry have been defined, including CD21/CR2, MHC class II and R9AP, much remains to be learned about host factors that support post-entry steps of EBV infection. To gain insight, we conducted a human genome-wide CRISPR-Cas9 screen. In addition to known EBV B cell receptors/co-receptors, we identified multiple new factors and pathways essential for EBV B cell infection. These included multiple B cell signaling pathways, actin cytoskeleton and nuclear import factors. Knockout of CD19, its tetraspanin chaperone CD81 or karyopherin subunit alpha 1 (KPNA1, also called importin-5) significantly impaired establishment of EBV infection. EBV and CD19 co-localized at early stages of EBV infection, and the CD19 cytoplasmic tail was important for EBV infection, suggesting CD19 signaling supports EBV uptake. Depletion of KPNA1 significantly impaired establishment of EBV infection, but did not impair production of infectious EBV upon reactivation of latently infected B cells. Together, these results provide insights into pathways that support EBV B cell infection and identify potential therapeutic targets.

Filoviruses, particularly those of the Orthoebolavirus genus, pose ongoing public health threats due to their increasing frequency and geographic spread. However, research has been impeded by the biosafety level 4 classification of filoviruses. We previously reported the generation of biologically contained Ebola, Marburg and Sudan virus as lower biosafety level-compatible filovirus systems. In the present study, we expanded this repertoire to include biologically contained Tai Forest, Bundibugyo and Reston virus, thereby creating a near-complete toolkit for currently recognized human-relevant orthoebolaviruses. VP30-deficient viruses were generated with matching VP30 expressing cell lines. More specifically, we developed and optimized a dual-reporter system in VeroE6, Huh-7 and A549 cells, combining a virus-encoded enhanced green fluorescent protein reporter as readout for viral replication with a stably cellular-expressed nuclear mCherry marker for cytotoxicity assessment. Next, we screened two repurposing-oriented compound libraries comprising 640 small molecules against Tai Forest, Bundibugyo and Sudan virus, and subsequently cross-validated against Reston and Ebola virus. This approach identified multiple candidates with broad-spectrum activity across orthoebolaviruses, while also revealing virus-specific antivirals, with robust activity observed in both primate- and human-derived cell lines. Together, this work establishes a versatile and experimentally tractable lower biosafety level-compatible platform for the study of human-relevant orthoebolaviruses and the systematic discovery of broad-spectrum antiviral countermeasures against filoviruses.

The rise of antibiotic-resistant bacteria demands alternative therapeutic strategies, with bacteriophage (phage) therapy and phage-derived enzybiotics emerging as promising approaches. However, identifying candidate phages against specific pathogens has historically been a bottleneck due to the need for cultivation methods to assess host range and lytic activity. Advances in metagenomic sequencing and the emergence of large-scale viral genome databases now provide an opportunity to accelerate this process computationally. Here, we present a large-scale mining of the MetaVR database to identify phages targeting human bacterial pathogens. By integrating direct host associations with CRISPR-spacer evidence, we linked 196,472 high-quality and complete viral genomes, representing 42,360 vOTUs, to 618 species of pathogenic and opportunistic bacteria. Functional enrichment analysis revealed distinct genomic signatures with viral lifestyle and host-range breadth: virulent phages were enriched in replication and structural functions, whereas temperate and broad host-range phages were enriched in anti-defense and regulatory modules. To characterize their lytic potential we annotated lysis-related protein families and their structural diversity, identifying 76 structurally novel lysis-associated proteins, including candidates targeting WHO priority pathogens. Focused analysis of endolysins revealed 592 structural clusters, with extensive sharing of endolysin repertoires among ESKAPE pathogens, suggesting candidates for broad-spectrum enzybiotic development. Selection analysis identified 167 endolysin families with sites under positive selection within functional domains, highlighting evolutionary diversification potentially associated with phage-host interactions. Together, our results establish a large-scale framework for connecting human bacterial pathogens to phages and their lytic machinery, providing a resource for prioritizing phage therapy and enzybiotic development.

Host genetic markers often fail to resolve regional origins in highly connected or panmictic marine species. The Japanese jack mackerel, Trachurus japonicus, is a commercially important fishery species around Japan that shows little or no detectable population structure. Here, we used nationwide multi-omics profiling to compare host genomic variation and gut microbiome composition in wild T. japonicus collected from coastal regions across Japan. We generated MIG-seq data for 43 individuals and 16S rRNA gene profiles for 24 individuals; after quality filtering, 19 individuals remained for matched host-microbiome comparison. Genome-wide host SNP analyses showed weak or absent geographic population structure, consistent with previous evidence of panmixia in Japanese waters. In contrast, gut microbiome composition showed geographic structuring based on Bray-Curtis dissimilarity and PERMANOVA, and this pattern was not explained by proximity to river mouths or host-related variables. Locality- or individual-associated bacterial lineages contributed to the observed differences in the microbiome, while chloroplast-associated and Cyanobacteria-assigned ASVs suggested recent dietary or environmental input. These results indicate that gut microbiome can show regional biological variation not apparent from host genetic markers alone. Our study provides a proof-of-concept example of integrating host genomics and gut microbiome profiling to evaluate regional characteristics and origins in highly connected marine animals.